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Thin Disk Models of Anomalous X-ray Pulsars 

K. Yavuz Ek§i'' and M.Ali Alpar"" 

""Sabanci University, Orhanli-Tuzla, Istanbul 34956, Turkey 

We discuss the options of the fall-back disk model of Anomalous X-Ray Pulsars (AXPs). We argue that the 
power-law index of the mass inflow rate during the propeller stage can be lower than those employed in earlier 
models. We take into account the effect of the super-critical mass inflow at the earliest stages on the inner radius 
of the disk and argue that the system starts as a propeller. Our results show that, assuming a fraction of the mass 
inflow is accreted onto the neutron star, the fall-back disk scenario can produce AXPs for acceptable parameters. 



1. Introduction 

The models for explaining the X-ray energy 
output of Anomalous X-ray Pulsars (AXPs, see 
PP for a review) can be classified in two groups. 
Magnetar models piSll] assign the persistent 
emission and bursts |5I6| to the extra-ordinary 
magnetic field of the neutron star while accretion 
models |7I8I9I10I11I12| employ neutron stars with 
intrinsic properties similar to other neutron stars 
and try to explain the X-ray emission and spin 
properties by the way matter is supplied. AXPs 
are very similar to Soft Gamma Ray Repeaters 
(SGRs) 13. 14 15. and probably have the same 
nature. 

Magnetar model is successful in modelling the 
SGR/AXP bursts and spin-down rates but can 
not explain the period clustering [16| except in 
one set of field decay models under special condi- 
tions |17| . Interaction of the neutron star with a 
fall-back disk towards equilibrium can naturally 
account for the period clustering. 

Mereghetti & Stella (1995) Y proposed that 
AXPs were accreting from a very low mass com- 
panion. There are severe upper limits for the 
mass of such a companion and no orbital motion 
signatures have been detected. Van Paradijs et al. 
(1995) |H| proposed a fossil disk model in which 
the disk results from the final spiral-in destruc- 
tion of the binary. Such a model addresses the 
absence of the binary companion but can not ac- 



count for the young ages of AXPs indicated by the 
supernova remnant (SNR) associations |18I19| . 
Fall-back disk models |9ll0lllir^ can address the 
young ages of AXPs and the absence of a mass 
donor companion. In all fall-back disk models 
the spin down of the neutron star is achieved by 
the propeller effect |^( 11:^11:^:^1:^31:^41^^ in which 
the rapid rotation of the neutron star prevents 
inflowing matter from accreting onto the surface 
of the star. Even when the star is slow enough 
that accretion can commence, it can continue to 
spin-down until torque equilibrium. AXPs must 
be those spinning-down accreting systems near 
torque equilibrium. 

Recent detections in the infrared (and optical) 
support the existence of a disk |27|. The magne- 
tar models do not have a prediction for infrared 
emission. Fall-back disks are not necessarily geo- 
metrically thin but, as the only models available 
for calculation, they warrant a careful study to 
decide whether interaction with a fall-back disk 
can produce neutron stars with AXP properties 
at the ages indicated by SNR-AXP associations. 
In this work we discuss the options of the fall-back 
disk model in view of the recent disk solution em- 
ployed by Ek§i & Alpar (2003) [E]. 

2. Fall-back Disk Models 

Alpar (2001) JOI proposed a classification of 
young neutron stars in terms of the absence or 
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presence and properties of a fall-back disk. Ac- 
cording to this model AXPs, SGRs and Dim 
Thermal Neutron Stars (DTNs) have similar peri- 
ods because they are in an asymptotic spin-down 
phase in interaction with a fall-back disk. The 
different classes represent alternative pathways of 
neutron stars. Radio pulsars have no disks or 
encounter very low mass inflow rates while Ra- 
dio Quiet Neutron Stars (RQNS) have such high 
mass inflow rates that their pulse periods arc ob- 
scured by the dense medium around. This model 
employs constant mass inflow rate corresponding 
to a constant equilibrium period. As a fall-back 
disk is not replenished, the mass inflow rate of the 
disk should decline while the disk diminishes. The 
mass flow rates characterizing different evolution- 
ary paths given in this model are then represen- 
tative values and are to be replaced by the initial 
mass of the disk in a time-dependent model. 

Chatterjee et al. (2000) (CHN ,9,) employed a 
time-dependent thin disk model |28I29| in which 
the angular momentum of the disk is constant and 
the mass of the disk declines with a power-law. 
Accordingly, the mass inflow rate also declines as 
a power-law 



-to 



(1) 



where Mq is the initial inflow rate and Tq is taken 
to be the dynamical time scale, ~ flj^^ at the 
initial inner radius of the disk. The power-law in- 
dex a depends on the opacity regime. For opaci- 
ties of the form k — kqp'^T^, thin disk equations 
can be solved to obtain 

_ 18a -Ab + 38 
""17^^26 + 32- ^' 

This yields, a = 19/16 for electron scattering 
opacity (a — b = 0), and a — 5/ A for bound- 
free opacity (a = 1 and b — —7/2). CHN as- 
sumed electron scattering to be dominating, but 
employed a = 7/6 for analytical convenience. 
As the mass inflow rate declines, the equilib- 
rium period increases and the system can never 
reach a true equilibrium. AXPs are those sys- 
tems accreting from a diminishing disk while spin- 
ning down towards an ever receding equilibrium 
(quasi-equilibrium stage) . 



Francischelli & Wijers (2002) ^ST criticized the 
CHN model, noting that whether the disk model 
can produce AXPs or not is sensitively depen- 
dent on the value of the power-law index a. 
They showed that after ~ 20 years, the disk is 
bound-free opacity dominated and for the cor- 
responding a, the model did not produce AXPs 
with the nominal disk mass employed by CHN 
(Mo = 0.006Mq). They also criticized the torque 
model employed in the CHN model, arguing that 
it is too efficient and showed that less efficient 
torque models can not produce AXPs (see also 

m)- 

While FracischeUi & Wijers (2002) [H] is right 
in arguing that a is critical in determining the 
period evolution for a single initial disk mass, it 
is not possible to reject the fall-back disk model 
on this ground. For every a, one can find an ap- 
propriate initial mass flow rate Mq (correspond- 
ing to an initial disk mass Mq) that leads to a 
period evolution in which the system enters the 
quasi-equilibrium stage at P = 5 s, producing 
an AXP. The time at which the transition to the 
quasi-equilibrium stage occurs, Ttr, and the mass 
flow rate at this moment, Mtr, remains nearly 
the same for each a and its corresponding Mq. 
One can obtain AXPs with the CHN model by 
employing a = 5/4 corresponding to bound-free 
opacity as long as one assumes an initial disk mass 
of Mq = O.O69A/0. This larger initial disk mass 
does not lead to higher luminosity at the accre- 
tion stage because the disk decays more rapidly 
with the corresponding a. 

Ek§i & Alpar (2003) (EA03 [E]) emplo yed an- 
other time-dependent thin disk solution j28l33| 
for the evolution of the disk during the propeller 
stage. In this solution mass of the disk is constant 
while the angular momentum of the disk increases 
due to the torque 

-/3 



J^Jo 



To 



at the inner boundary of the disk where 
^ 14a + 22 



15a -26 4- 28 



(3) 



(4) 



For electron scattering opacity f3 — 11/14 and for 
bound-free opacity dominated disks (3 = 18/25. 
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In this solution mass flows outwards, but Pringle 
(1991) 33 numerically shows solutions in which 
the mass flows inwards, and a torque as given in 
Eq. ||2Jl is necessary to stop accretion at the in- 
ner boundary. Pringle (1981) (34| mentioned that 
"such a solution might represent a disc around 
a magnetized star which is rotating sufficiently 
rapidly that its angular velocity exceeds the Ke- 
plerian velocity at the magnetosphere" - i.e. a 
propeller. Assuming the inner radius Rm scales 
with the Alfven radius Ra oc Af ^^Z*" and referring 
J oc ^JGMRmM (X Af^/*", one can calculate an 
effective a for the propeller regime. For electron- 
scattering opacity one obtains ap = 11/12 and 
for bound-free opacity ap = 21/25. EA03 used 
another scaling for determining the inner radius 
in the propeller regime Rm oc M~^/^ J35;. This 
implementation for the inner radius yields J oc 
]\^9/io g^j-^^ Qj^^ obtains ap = 55/63 for the elec- 
tron scattering opacity and ap = 4/5 for the 
bound-free opacity. 

For the sake of consistency with the constant 
disk mass solution, EA03 assumed that the in- 
flowing matter flung away by the magnetosphere 
of the neutron star returns back to the disk some- 
where away from the inner boundary of the disk. 
The mass flow rate in the propeller type of solu- 
tion declines with a softer power-law index than 
it does in the accretion type of solution because 
in the former case matter cannot accrete and the 
decrease in the flow rate is only by the viscous 
spreading of the disk. Conservation of the mass of 
the disk during the propeller stage makes it possi- 
ble to obtain AXPs with a 50 times smaller initial 
disk masses than required by the CHN model. 
Again, this does not mean that mass flow rate 
during the accretion stage will be smaller by this 
order of magnitude. 

The luminosities of AXPs being less than Lx — 
7 X 10^^ erg s~^ imply that less than Mx — 
4 X 10^^ g s~^ of matter is accreting onto these 
objects. As Menou et al. (2001) (3^1 have shown, 
fall-back disks with M < 10^^ g s~^ are neu- 
tralized and can not be accreting. This means 
that AXPs must have mass flow rates in excess of 
~ 10^^ g and only a fraction r] of the inflow- 
ing matter should accrete onto the neutron star. 
In CHN model Mtr ^ 10^^ g and one needs 



to assume that a fraction 77 = 0.1 of the inflowing 
mass accretes onto the neutron star, in order to 
explain the observed luminosities of Lx S 10'^^ 
erg s~^. In EA03 model one needs to assume 
rj ~ 0.01 in order to explain the luminosities be- 
cause the system reaches the quasi-equilibrium 
stage rather rapidly, at an epoch during which 
the mass flow rate is still high. The difference 
is due to the different formulations employed for 
the magnetic radius of the disk during the pro- 
peller stage which in turn effects the magnitude 
of the torque and the time the system reaches the 
quasi-equilibrium stage. 

3. The Model Equations 

For finding the inner radius, we follow the 
EA03 model The system initially is in the 

super-critical propeller regime and the inner ra- 
dius of the disk is determined by 

Rm = {fi^Kes/8TTV2cn^y/^ 

= 31i?*Ai3^i^(Po/15ms)i/s (5) 

while the corotation radius is Rc — 
10i?*(P/15ms)^/'^. The dynamical time scale 
Td = 1/flKiRm) correspondingly is Td ~ 
1.23 X 10-V3o'^(^o/15ms)i/4 which is one or- 
der of magnitude greater than that employed 
by CHN. Greater Tg makes it possible to ob- 
tain the same mass flow history (and same Mtr) 
starting with a smaller Mq by (and smaller Mq). 
The calculations of Pringle (1991) 33 show that 
the power-law evolution starts after a viscous 
time scale rather than the dynamical time-scale 
employed in early models. Considering the am- 
biguity in To, one should not attempt to be very 
definite on the initial mass of the disk in models 
employing the power-law evolution. 
We model the disk torque as 

k - VGMRmM{l - uj,/uj,g) (6) 

where ujeq is the equilibrium fastness parameter 
which we take to be 0.85 in this work. As noted 
by 12] and [22], this is a very efficient propeller 
torque model compared to some other available 
propeller torques obtained through simple energy 
or angular momentum arguments 21 22 23 . Re- 
cently Ikhsanov |26| revived the propeller model 
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of Davies & Pringle (1981) 24 which incorpo- 
(X — oj*. This model agrees with Eq. © 
in the large limit. In the form we adopt, the 
torque entails asymptotic approach to rotational 
equilibrium which is necessary for addressing the 
period clustering of AXPs. The recent numeri- 
cal work [HKj also indicate an efhcient propeller 

torque (j* a — ri^"^)- 

Being aware that modifications in a can be 
compensated by modifying the initial inflow rate 
which in turn depends on the uncertain value of 
the dynamical time scale, we simplify the EA03 
model as follows: 

a = i'/' '"''^^'-^ (7) 
[4/5 if > 1. ^ ' 

The value 5/4 is motivated by the accretor type 
of solution and the value 4/5 is motivated by 
the propeller type of solution. Both values corre- 
spond to the bound free opacity regime. The elec- 
tron scattering opacity prevails only in the 20 
years and does not effect the period evolution sig- 
nificantly. 

4. Results and Discussion 

We have followed the evolution of the star-disk 
system for 10^ years, solving the torque equa- 
tion lO numerically with the adaptive step-size 
Runge-Kutta method [S7] . 

In Figure ^ we show five evolutionary tracks 
on the P — P diagram that can cover seven 
AXP/SGRs for which period derivatives are avail- 
able. Figure (0) shows the corresponding pe- 
riod evolutions for each of the tracks in Fig- 
ure 1^, tracks with the same numbers having 
the same parameters. We assumed the radius 
and mass of the neutron star to be i?* = 10® 
cm and Af, — 1.4M0, respectively. Initial period 
is Pq = 15 ms and we have assumed that the disk 
neutralizes when the mass flow rate drops below 
IQi® g s-i 

Extending the analysis to disk evolution in 
the propeller regime provides favorable effective 
power-law indices < 1. However, in the avail- 
ability of choosing a larger initial disk mass one 
can always obtain AXPs even if the disk looses 
mass with a = b/A. 



Considering the effect of the radiation pressure 
on the inner radius indicates a larger dynamical 
time scale, Tq. With a larger Tq, one can produce 
the same mass flow history with a smaller initial 
mass inflow rate corresponding to a smaller ini- 
tial disk mass. The numerical calculations (see 
|33j) indicate that the power-law evolution cor- 
responding to the self-similar solutions start in a 
viscous time-scale which is still larger than the 
dynamical time scale. Given the ambiguity in Tq, 
estimations of initial disk mass from power-law 
models discussed above, should not be taken very 
seriously. 

Modifications in a and Tq can always be com- 
pensated in the sense that one can obtain the 
same period evolution by tuning the initial mass 
fiow rate. Such modifications then lead to the 
same mass fiow rate at the quasi-equilibrium 
stage at which the system is observed as an AXP. 

The part of the model that is critical in deter- 
mining whether fall-back disks can produce AXPs 
or not is the propeller torque as first shown by 
|31j . If the torque model is inefficient, one must 
then employ a high mass fiow history to spin- 
down the neutron star to AXP periods which in 
turn results with very high luminosities in the 
quasi-equilibrium stage. As the torque is depen- 
dent on the inner radius of the disk, the results are 
effected by the uncertain factors (^ = Rm/RA) in 
determining the inner radius. 

It is possible that a fall-back disk loses mass 
also in the propeller regime because of the pro- 
peller effect itself, depending on the fastness pa- 
rameter. The simple power-law models for the 
fiow rate are working models which do not take 
into account the coupling of the star to the struc- 
ture of the disk. A more realistic treatment of fall- 
back disks would consider the effects of the mag- 
netic field on disk evolution, realistic opacities, 
including the iron line opacity |SH| for the fall- 
back disk which is rich in heavy elements, and the 
dependence of the accretion rate itself on the mo- 
tion of the inner radius M = 2'KRmX'{Rm)[Rm — 
Vr{Rm)] (see ISni)- All of these effects would re- 
sult in time evolution of the disk that cannot be 
described by simple power laws. Taking account 
of the complexities of the real problem is not 
likely to produce qualitative disagreement with 
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the evolutionary timescales as long as the torque 
is efficient. 
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Figure 1. Evolutionary tracks in the P — P di- 
agram covering all AXPs and two of the SGRs 
with known period derivatives. The thin solid 
lines show the rapid spindown phase, the thick 
solid lines stand for the "bright AXP" stage and 
the dashed lines show the "dim AXP" stage after 
the disk is neutrahzed (M < 10^'' g s~^). The 
dots denote the radio pulsars. AXPs are shown 
by filled circles and SGRs by squares. With dif- 
ferent /i3o and Mq, the model can cover all AXPs 
and SGRs. The parameters for each track are 
as follows: (l)Pi2 = 1.8, Mq = 5 x lO^^ g s'^; 
(2)Pi2 = 2.8, Mo = 1 X 1028 g s-i; (3)Pi2 = 5.5, 
Mo = 5 X 10^^ g s-i; (4)Pi2 = 8, Mo = 3.2 x lO^^ 
g s-i; (5)Pi2 = 6, Mo = 2 X 10^^ g The ini- 
tial period is 15 ms for all simulations. 
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Figure 2. Period evolution of AXPs. The thin 
solid lines show the rapid spindown phase, the 
thick solid lines stand for the "bright AXP" stage 
and dashed lines show the "dim AXP" stage af- 
ter the disk is neutralized. The dashed and dotted 
horizontal lines show the period range of AXPs 
and SGRs. The parameters are the same as in 
Figure and the numbers labelling the tracks 
correspond to the numbers of the tracks in Fig- 
ure (QJ. 



